Fermi gamma-ray 'bubbles' from stochastic acceleration of electrons 
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Gamma-ray data from Fermi-LAT reveal a bi-lobular structure extending up to ~ 50° above and 
below the galactic centre, which presumably originated in some form of energy release there less 
than a few million years ago. It has been argued that the 7-rays arise from hadronic interactions 
of high energy cosmic rays which are adverted out by a strong wind, or from inverse-Compton 
scattering of relativistic electrons accelerated at plasma shocks present in the bubbles. We explore 
the alternative possibility that the relativistic electrons are undergoing stochastic 2"'*-order Fermi 
acceleration by plasma wave turbulence through the entire volume of the bubbles. The observed 
7-ray spectral shape is then explained naturally by the resulting hard electron spectrum modulated 
by inverse-Compton energy losses. Rather than a constant volume emissivity as in other models, we 
predict a nearly constant surface brightness, and reproduce the observed sharp edges of the bubbles. 

PACS numbers: 98.70.Rz, 98.35.Jk, 95.30.Qd, 98.70.Sa 



Recent data from the Fermi-LAT satellite has re- 
vealed [T1I2] the presence of giant 7-ray lobes ~ 40° wide, 
extending up to ~ 50° above and below the galactic cen- 
tre (GC). The energy spectrum of the emission from these 
'Fermi bubbles' is dN/dE ^E'"^ from - 1 - 200 GeV, i.e. 
considerably harder than conventional foregrounds. Fur- 
thermore, the bubbles exhibit an almost constant surface 
brightness with hard edges. While the template subtrac- 
tion technique used to reveal the bubbles may not be ap- 
propriate at these high energies, the resulting systematic 
effects are not easy to assess. However the bubbles do cor- 
relate with features at other wavelengths, viz. data from 
the ROSAT X-ray satellite [3] show a limb-brightened, 
conical structure close to the galactic plane which coin- 
cides with the edges of the Fermi bubbles. The bubbles 
also line up with a claimed excess in microwaves at lower 
galactic latitudes — the so-called 'WMAP haze' [J]. 

Although extended lobes have long been seen in other 
galaxies in radio. X-rays and 7-rays, their presence in the 
Milky Way is surprising. There is no radio emission from 
these bubbles, unlike those seen in the majority of active 
galaxies. Moreover their morphology (symmetry with re- 
spect to the galactic plane and alignment with the GC) 
suggests that the central supcrmassive black hole is the 
energy source. However it is supposedly in a quiescent 
state so it is a puzzle how the bubbles have formed; un- 
derstanding this would provide an excellent probe of this 
region which is otherwise obscured by the galactic disk. 
The bubbles may play an important role in the dynamics 
of our galaxy and constitute a source of cosmic rays (CR). 
While they are prominent at high galactic latitudes, the 
associated signal close to the plane, while uncertain, con- 
stitutes a background for indirect dark matter searches. 
It is therefore important to understand and model the 
origin of the non-thermal emission from the bubbles. 

While the mechanism responsible for the formation of 
the bubbles is not necessarily the same as the source of 
the 7-ray emission today, it is useful to recall their gen- 
eral properties. The limb-brightened shell in the ROSAT 
data implies a shock front at the bubble edges, but from 



the observed cavity hot low density gas is inferred to 
fill the bubble interiors. Assuming a low density {n ^ 
10~^cm~'^) gas at T 2keV and shock velocities U < 
1000 km s^^, the energy is estimated to be ~ lO^''^"^^ erg 
in hot gas and the age to be ~10''(?7/1000kms~^)yr [1]. 
Suggested mechanisms for providing such an energy on 
this timescale include jets emanating from the central 
black hole [5], star forming regions close to the GC [6] or 
repeated star accretion onto the central black hole 7]. 

The observed 7-rays may be generated by hadronic in- 
teractions of high-energy CR protons or nuclei {i.e. 
decay) provided that the ambient gas-density is not too 
low. It has been proposed |6] that protons and nuclei ac- 
celerated by supernova remnants (SNRs) in star-forming 
regions very close to the GC could be adverted by a 
strong wind out to kiloparsec distances above the plane. 
If the confinement time is larger than all other timescales, 
the hard power law spectrum of the 7-rays would simply 
reflect the source spectrum of the protons. The spectral 
shoulder at 1 GeV can be explained by the pion bump. 

Another possible mechanism is the inverse-Compton 
(IC) scattering of high energy electrons off ambient radi- 
ation fields (CMB, far infra-red (FIR) and optical/UV). 
The spectral feature seen at a few hundred GeV may 
reflect a cut-off in the electron spectrum at a similar en- 
ergy, either due to energy losses or due to the competition 
between an energy-dependent acceleration rate and the 
finite age of the bubbles. Furthermore, the WMAP haze 
may well arise from synchrotron radiation of these elec- 
trons in the ambient magnetic field. A crucial question 
then is how are the electrons accelerated. 

The standard paradigm for the acceleration of galactic 
CRs is diffusive shock acceleration (DSA) by the l^*-order 
Fermi process, which predicts power-law source spectra 
with index close to —2. There are at least four regions 
where shocks may be present: at the GC, inside a jet 
emanating from the GC, at its termination shock at the 
upper/lower edges of the bubbles, and at the shocked ex- 
terior of the bubbles. So far there is only evidence from 
ROSAT data for a shock at the bubble exterior. In any 
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case, presuming difFusive-convective transport from the 
acceleration site through the bubble volume, it is diffi- 
cult to see how the electrons can maintain their hard 
source spectrum. The energy loss time due to IC scat- 
tering for the ©(TeV) energy electrons present through- 
out the bubble is only a few times 10^ yr; however even 
with a convection velocity as high as u ~ 1000 km s~^, it 
would take the electrons 10^ yr to cross the required dis- 
tance of 0(10) kpc. The leptonic source model [7] there- 
fore invokes hundreds of consecutive shocks in order to fill 
the whole bubble with freshly accelerated electrons. This 
would however imply a constant volume emissivity which 
in projection would yield a characteristic bump-likc pro- 
file with soft edges, in contrast to what is observed. 

We consider instead the stochastic acceleration of high 
energy electrons by isotropic, large-scale turbulence in 
magnetosonic waves [8J. Such 2"'^-order Fermi accelera- 
tion accounts well for the radio emission from supernova 
remnants [HI HO] and the extended lobes of radio galax- 
ies [11] , and may even be the acceleration mechanism for 
ultra- high energy cosmic rays [12]. The shock front at 
the bubble edges suggests that they may have been pow- 
ered by a jet emanating from the massive black hole at 
the GC that was active a few million years ago. MHD 
modelling [5] of a two-component plasma explains the for- 
mation of a bubble by a light but over-pressured jet with 
~ 16 % of the Eddington luminosity, and also predicts a 
shock coincident with the ROSAT shell. Plasma instabil- 
ities, in particular Rayleigh- Taylor and Kelvin-Helmholtz 
instabilities, would then generate turbulence at the outer 
shock that is convected into the bubble interior by the 
downstream plasma flow. The free energy dissipation 
rate Q = Cipu^/L is determined by the scale of turbu- 
lence injection, L, and the eddy velocity at the injection 
scale, u = Vedd{L), where Ci = 0.485 is the 1-dimensional 
Kolmogorov constant [10] . The energy density at scale k 
is then given by W{k) = (u^/47r)L~^/^fc~^^/^. Applying 
the Rankine-Hugoniot conditions at the shock, the eddy 
velocity at the injection scale, u, and the magnetosonic 
phase velocity, vf, vary with the distance x = (,L from 
the shock as [TU] : 
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the dissipation scale to be the mean-free path, thus ren- 
dering the spatial diffusion coefficient D^x = ldc/3 en- 
ergy independent. If additional small-scale turbulence 
is present (possibly responsible for spatial diffusion [S^), 
then the mean-free path can be smaller. 

The temporal evolution of n{t,p) dp, the number den- 
sity of electrons with momentum between p and (p + dp), 
is dictated by the Fokker-Planck equation 13j, 
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where the diffusion coefficient in momentum for scatter- 
ing by fast magnetosonic waves is [B]: 
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The second term in Eq. |3] describes diffusion in momen- 
tum as well as systematic energy gains on the charac- 
teristic timescale ^acc ~ P^/Dpp, which is also energy 
independent. Diffusive losses from the acceleration re- 
gion can be accounted for by escape on the timescale 
tcsc — /Dxx- Finally, the electrons lose energy through 
IC scattering and synchrotron radiation which are both 
accounted for by the energy dependent cooling time 
^cooi = -p/{dp/<it). 
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where U is the shock velocity, va the Alfven velocity 
(which we assume to be constant and equal to the speed 
of sound Vs.o at the shock) and a = 3 — 16v^q/U^. 

At small enough scales Id = 1/fcd = L{va/u)^, the ki- 
netic energy of the turbulence becomes comparable to the 
magnetic field energy, Wcdd(^d) ~ resulting in transit- 
time damping. With parameters to be justified below, it 
turns out that for all energies of interest, the gyro-radius 
Tg of the electrons is always smaller than this dissipation 
scale Id such that gyro-resonant interactions with magne- 
tosonic turbulence are not possible. Therefore, we adopt 



FIG. 1. Relevant timescales (top) and the electron spectrum 
(bottom), at various distances x = from the shock. 

Because of the energy-independent spatial diffusion co- 
efficient, the so-called "hard-sphere" approximation [14] 
is exact which makes the problem amenable to analyti- 
cal solution. If the escape rate is not much bigger than 
the acceleration rate, i.e. tacc ^ ^esc, the steady state 
spectrum n{p) at a fixed position can be described as 
a power law with a spectral cut-off above (and pile- 
up around) a characteristic momentum poq, defined by 
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The spectral index, — cr = 1/2 — ^9/4 + iacc/^csc, is de- 
terniined by the ratio of acceleration and escape times, 
and asymptotically approaches —1 as iacc/iesc 0. 

Anticipating that the acceleration time is smaller than 
the lifetime of the bubbles, we justify the use of the 
steady-state solution for acceleration volumes that are 
being advected with the downstream plasma. It is suf- 
ficient to consider the variation of the acceleration and 
escape times with the distance from the shock, which de- 
termines the spatial dependence of the electron spectrum. 
This hierarchy of timescales assures that the variation 
with position happens adiabatically, such that the elec- 
trons can always relax to their steady-state spectrum. In 
the upper panel of Fig.[l]we show the different timescales 
in the problem as a function of energy for the parameters 
discussed below. Although icooi is of the same order as 
the dynamical time tconv around 10 GeV, we expect that 
the steady state spectrum is reached in a time t ^ iacc, 
as has been shown explicitly [TB] for ionisation losses. 

The relative normalisation of the electron spectrum is 
fixed by noting that the total energy in relativistic elec- 
trons at any position is a constant fraction of the free 
energy dissipated along with the downstream plasma up 
until this position. This does not however fix the abso- 
lute normalisation which depends on the microphysics of 
the acceleration process, in particular the injection mech- 
anism. We determine the 7-ray volume emissivity due to 
IC scattering off the CMB, FIR and optical/UV back- 
grounds adopting the interstellar radiation fields from 
GALPROP at a reference height of 4kpc above the GC. 
For the parameters discussed below we show the electron 
spectrum E^n^- for different distances from the shock in 
the lower panel of Fig. [T] 

We now discuss the parameters that can repro- 
duce the observed 7-ray flux — both its spectrum 
and morphology. Kelvin-Helmholtz instabilities have 
been observed to be generated on kpc scales in MHD 
simulations of the Fermi bubble gas [S], so we choose 
the scale of turbulence generation to be L = 2 kpc. 
The shock velocity can in principle be determined 
kinematically from the variation of its position with 
time (the shock needs ^ 50 (C//10^ cms~^) yr to move 
a distance corresponding to the 1" resolution of the 
Chandra X-ray observatory) or possibly inferred from 
the observed shock heating. We fix [/ = 2.6 x 10^ cms^^, 
a value consistent with MHD simulations [5]. Finally 
the Alfven velocity is given by the square root of 
the ratio of magnetic field energy density to thermal 
plasma energy density: /3a — va/c — yJUs/Up. Hence 
/3a > 2.8 X 10"'' for an estimated upper limit on the 
thermal gas density n < 10~^ cm~^ [1] and a magnetic 
field B = 4^G (suggested by radio observations of 
the edge-on spiral galaxy NGC 891 [H]). We adopt 
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FIG. 2. The average 7-ray flux from the Fermi bubbles [T] 
compared to the spectrum from our model (the contributions 
from inverse-Compton scattering on the CMB, FIR and op- 
tical/UV backgrounds are shown from left to right as dot- 
dashed lines). The spectra from a hadronic ^ (dashed line) 
and a leptonic DSA model ii?^ (dotted line) are also shown. 



/3a = 5 X 10 The gyro radius of relativistic elec- 
trons is then - 7.5 x {B / 4, ^iG)-^ {E / GeV) cva., 
which is much smaller than the dissipation length 

> 8 X 10i9(L/kpc)(C//10«cms-i)-3(/3A/10-3)3cm 
even for 0(10) TeV electrons, thus confirming the 
energy-independence of the acceleration and escape 
time. With these parameters we find a total energy in 
electrons above 100 MeV of ~ 10^^ erg which is over five 
orders of magnitude smaller than the required energy in 
protons in the hadronic emission model [B]. 

In Fig. [2] we show our predicted flux E'^J^ of high 
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FIG. 3. The 7-ray intensity as a function of distance from 
the bubble edge [1] is compared with our model predictions 
at 2 GeV (soUd), 10 GeV (long dashed), and 500 GeV (dot- 
dashed). The dotted line indicates the expected profile for 
both the hadronic model [B] and the leptonic DSA model [7]. 
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FIG. 4. Our model prediction for the synchrotron flux at 
{i, b) = (0°, 25°), for two assumed magnetic field values, com- 
pared to the inferred spectrum of the WMAP haze [4|. 



energy 7-rays (averaged over the surface of the Fermi 
bubbles) as a function of energy, and compare it to the 
data [T] as well as to the hadronic [3] and leptonic DSA [7] 
models. Note that our hard electron spectrum nicely re- 
produces the spectral shoulder around a GeV and the 
cut-off at a few hundred GeV. We show how the total 
gamma-ray flux is made up of contributions from IC scat- 
tering on the CMB, FIR and optical/UV backgrounds 
(dot-dashed line, from left to right). Since we expect the 
FIR and optical/UV contributions to decrease rapidly 
with distance from the disk, the emission from high lati- 
tudes should cut off above tens of GeV — a potential test 
of the model. In Fig.|3]we compare the data with the pre- 
dicted 7-ray intensity at 2 and 10 GeV as a function of 
distance from the bubble edge (calculated as in Ref. [T], 
i. e. averaging over arcs of great circles converging at the 



bubble centre). This test has not been done before; it 
is seen that our model matches the almost constant in- 
tensity in the interior with hard edges, in contrast to the 
hadronic [B^ and leptonic DSA [7 models. We also show 
our model prediction for higher (500 GeV) 7-ray energies; 
the presence of sharp edges in 7-rays can be tested by the 
forthcoming Cherenkov Telescope Array. 

The energy dependence of the profiles reflects the spa- 
tial variation of the electron spectrum with distance from 
the shock (see also Fig. [ij . Close to the shock, the accel- 
eration time is relatively small such that the spectrum is 
very hard, Pcq is large, and the spectral bump and cut-off 
appear at high energies. Further away from the shock, 
tacc becomes larger and the spectrum softer, while the 
bump and cut-off move to lower energies. The emission 
of the highest energy 7-rays (due to the highest energy 
electrons) is thus localised close to the shock and results 
in the limb-brightening above a few hundred GeV. Inter- 
mediate energy 7-rays can be produced from both high 
and intermediate energy electrons which have a more ex- 
tended distribution, leading to a flatter intensity profile. 

While the 'WMAP haze' |1] has not been observed in 
polarised emission |19j and may just be an artefact of 
the template subtraction [20], it has been proposed as a 
physical counterpart of the Fermi bubbles p]. However 
as seen in Fig. |4j the expected synchrotron flux in our 
model is of the required amplitude only if the magnetic 
field is as strong as 15 //G, several kpc from the plane. 

The hadronic model predicts a detectable fiux of neu- 
trinos for the proposed Mediterranean km'^ neutrino tele- 
scope [6] . However the observed bubble profile disfavours 
this model (as well as the leptonic DSA model) and in- 
stead favours 2"'^-order Fermi acceleration of electrons, 
which would not generate any neutrinos. 
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